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TECHNICAL PUBLICATION 


ON THE RELATION BETWEEN SUNSPOT AREA AND SUNSPOT NUMBER 

1. INTRODUCTION 


Prior to 1874, the primary measure of solar activity was sunspot number, 1-4 in particular, Wolf’s 
relative sunspot number R, given as R = k (lOg +f), where g is the number of sunspot groups, /is the num- 
ber of individual spots, and k is a correction factor, which Wolf originally assigned the value of 1. In 1882, 
Wolf’s successors at the Swiss Lederal Observatory in Zurich changed somewhat the counting method for 
R and assigned a new value of 0.60 for k, which reduced the new observations to the old scale. The Swiss 
Lederal Observatory continued to provide Wolf or Zurich sunspot numbers through 1980. (Beginning 
in 1981 and continuing through the present, the international sunspot number is provided by the Royal 
Observatory of Belgium.) 

Because the definition of R is rather arbitrary, dependent upon accurate measurements of g and 
f which can vary considerably even with identical instruments and seeing conditions, a more objec- 
tive measure of solar activity was desired— sunspot area. Beginning in May 1874, the Royal Greenwich 
Observatory (RGO) began cataloging daily sunspot positions and areas, based on systematic observa- 
tions at Greenwich, England; Cape Town, South Africa; and Kodaikanal, India. 5 The data were obtained 
directly from photographs by means of reticules divided into small squares and measured relative to Sun 
center. Typically, the position of each sunspot was determined to 0. 1 degrees (latitude and central meridian 
distance), and the area of each spot was determined to 0.1 millionth of the visible solar hemisphere. The 
RGO dataset ended in 1976. 

To continue the sunspot positional/area record beyond 1976, United States Air Lorce/National Oce- 
anic and Atmospheric Administration (USAL/NOAA) Solar Optical Observing Network (SOON) observa- 
tions were employed. 6 Unlike the RGO photographic observations, the USAL/NOAA SOON observations 
are visually determined. A worldwide network of solar observatories (including Boulder, Holloman, Lear- 
month, Palehua, Ramey, San Vito, and occasionally, Culgoora— located so that 24-hr synoptic coverage 
can be achieved) comprises the SOON system. At each observatory, daily sunspot drawings on an 18-cm 
diameter image of the Sun are routinely made, with the scaling of positions and areas of sunspots being 
performed by hand using Stonyhurst disc overlays for positions and circle/ellipse overlays for sunspot 
area starting in 1981— grids were used to estimate sunspot area prior to 1981. Typically, positions are 
determined to 1 -degree accuracy and areas are determined by rounding to the nearest 10 millionths of the 
visible solar hemisphere. The complete set of observations from May 1874 through the present, both RGO 
and SOON, is available online at http://science.nasa.gov/ssl/PAD/SOLAR/greenwch.htm. 


The relation between total corrected sunspot area A (in millionths of the visible solar hemisphere) 
and sunspot number R is not particularly good for daily values, although it has often been stated that the 
relationship becomes better when using monthly or yearly averages. In fact, the claim 5 ’ 7-11 has been made 
that for monthly or yearly averages the average relation between A and R is simply described as A= 16.7 R. 

In this study, yearly averages of A, R, and AIR for 1875-2004 are examined to determine the valid- 
ity of the oft-quoted relationship between sunspot area and sunspot number. Also, because of the changes 
that have occurred in the solar activity record, the interval of 1958-1998 is examined in greater detail by 
comparing RGO and SOON sunspot areas against Rome Observatory sunspot areas and comparing Swiss 
Federal Observatory and Royal Observatory of Belgium sunspot number values against Rome Observa- 
tory determinations of sunspot number. Previously, it has been shown that both sunspot area and sunspot 
number, as determined today, appear slightly different (underestimated and overestimated, respectively) 
from earlier years when sunspot areas were determined by the RGO and sunspot numbers were deter- 
mined by the Swiss Federal Observatory. 12 
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2. RESULTS AND DISCUSSION 


2.1 Sunspot Areas, Sunspot Numbers, and Their Ratios (1875-2004) 

Figure 1 displays yearly plots of R, A, and AIR for the interval 1875-2004. The numbers at the bot- 
tom (12, 13, etc.) refer to specific sunspot cycles. The thick and thin vertical lines refer, respectively, to the 
occurrences of sunspot number maximums and minimums. The numbered triangles across the top (1,2, 
and 3) refer to times when changes occurred in the solar activity record, where No. 1 refers to the change 
in the method for counting sunspot number (using A: = 0.60) that occurred in 1882, No. 2 refers to the end 
of the RGO record and the beginning of the SOON record in 1976/1977, and No. 3 refers to the end of 
sunspot number determination at the Swiss Federal Observatory and the start of sunspot number determi- 
nation at the Royal Observatory in Belgium, specifically, the Solar Influences Data analysis Center (SIDC) 
in 1980/1981. The horizontal lines are the medians for each parameter. Also shown is the mean, standard 
deviation, and the results of a runs test for each parameter, where n a refers to the number of yearly values 
above the median, n b refers to the number of yearly values below the median, n r refers to the number of 
mns in the sample, and z is the statistic for determining randomness (where z > 1.96 indicates that the fluc- 
tuations are nonrandom in nature, in this case, being due to the systematic variation of the solar cycle). 13 

As can be seen, based on the yearly averages, A and R peak together for all cycles, except cycles 
20, 21, and 23, when A peaked 2 yr after R. For AIR, it was at minimum value at sunspot number mini- 
mum for cycles 13, 15, 16, 19, 20, and 22. It occurred 2 yr before sunspot number minimum for cycle 17; 
1 yr before sunspot number minimum for cycles 14 and 23; and 1 yr after sunspot number minimum for 
cycles 12, 18, and 21. AIR was at maximum value 2 yr before sunspot number maximum for cycle 16; 1 yr 
before sunspot number maximum for cycles 18 and 20; at sunspot number maximum for cycles 12 and 14; 
1 yr after sunspot number maximum for cycle 17; 2 yr after sunspot number maximum for cycles 19, 22, 
and 23; 3 yr after sunspot number maximum for cycle 21; 4 yr after sunspot number maximum for cycle 
13; and 5 yr after sunspot number maximum for cycle 15. 

For the entire interval, 1875-2004, AIR has a median value of 13.6 and a mean value of 13.1, with 
a standard deviation of 3.5. During the RGO timeframe (1875-1976) AIR averages 14.1, with a standard 
deviation of 3.2; and during the SOON timeframe (1977-2004) AIR averages 9.8, with a standard devia- 
tion of 2.1). Presuming that AIR is normally distributed during each of the subtimeframes, one finds that 
the t statistic for independent samples 14 is statistically meaningful, meaning that AIR is definitely lower 
during the latter timeframe than the former timeframe. 

Figure 2 depicts scatterplots of yearly values of A versus R for the two separate timeframes, 1 875— 
1976 (bottom) and 1977-2004 (top). Clearly, during each timeframe A is strongly correlated against R, 
although the slopes and y-intercepts are different. As an example, given a value of 100 for R, during the 
former timeframe it would suggest a value of A of =1 ,538 ± 174 (1-a accuracy), but only a value of A of 
=1,076 ± 123 during the latter timeframe. This suggests that yearly area averages are =30 percent smaller 
during the present timeframe than during the RGO timeframe. 
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Figure 1. Variation of yearly averages of sunspot number R (lower panel), sunspot area A 
(middle panel), and ratio AIR (upper panel) for 1875-2004. 




Figure 2. Scatterplots of A versus R for 1875-1976 (lower panel) and 1977-2004 
(upper panel). See Nomenclature list for details. 
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Figure 3 compares cycle 23 yearly values of R (bottom), A (middle), and AIR (top) against the 
mean parametric values of cycles 12-22 relative to the epoch of sunspot minimum year. Values of R for 
cycle 23 lie above the mean for all years, except yr 3, with peak value occurring in yr 4, as found for the 
mean of cycles 12-22. In terms of maximum amplitude, cycle 23 ranks fifth below cycles 19, 22, 21, and 
18, and above cycles 12-17 and 20. The behavior of cycle 23 ’s R is such that minimum for cycle 24 is 
anticipated for year 10, that is, 2006. 

Concerning cycle 23 ’s values of A, they lie below the mean of cycles 12-22 for all years, except 
the last three, when they are essentially the same as the mean of cycles 12-22, having a maximum yearly 
area that peaked in yr 6 that is only slightly larger than values for cycles 12 and 14. Concerning the ratio 
AIR, cycle 23 ’s ratio has always fallen below the mean of cycles 12-22, and, interestingly, its maximum 
ratio value is similar to what was seen for cycles 21 and 22, suggesting that something clearly is amiss 
with sunspot areas measured today as compared to years prior to 1977. 

2.2 Comparisons Against Rome Observatory Measures (1958-1998) 

Figure 4 (top panel) shows the number of days during each year for 1958-1998 that the Rome 
Observatory reported measurements of sunspot area, number of sunspot groups and number of individual 
spots, and the latter two parameters allowing for the computation of sunspot number as determined by the 
Rome Observatory. The bottom panel of figure 4 compares sunspot number s x and s R , where s x refers to 
the Swiss Federal Observatory (Zurich) timeframe for 1958-1980 and to the Royal Observatory of Bel- 
gium (SIDC) timeframe for 1981-1998 and s R refers to the Rome Observatory for the years 1958-1998, 
all values computed using the same days. The bottom-middle panel of figure 4 compares sunspot area 
and A r , where A x refers to the RGO timeframe for 1958-1976 and to the USAF/NOAA SOON timeframe 
for 1977-1998 and A R refers to the Rome Observatory for 1958-1998, again all values computed using the 
same days. The top-middle panel of figure 4 compares the ratios A x /s x and A R /s R for 1958-1998. 

Notice that s x and s R are virtually identical for 1958-1980, but differ markedly during the latter 
interval (1981-1998). For 1958-1980, the average deviation between Zurich sunspot number and Rome 
Observatory sunspot number measured only 3.7 units of sunspot number, with a standard deviation of 3 
units of sunspot number. For 1981-1998, the average deviation between SIDC sunspot number and Rome 
Observatory sunspot number measured 10.9 units of sunspot number, with a standard deviation of 12.5 
units of sunspot number. In particular, SIDC sunspot number for maximum amplitude for cycle 22 (in 
1989) measured 160.2, some 34.6 units higher than measured by the Rome Observatory sunspot number. 

Concerning sunspot areas: whereas Rome Observatory areas tended to be slightly below areas 
measured by the RGO (average deviation equal to 111.1 millionths of the visible solar hemisphere, with 
a standard deviation of 127.7 millionths of the visible solar hemisphere), for the USAF/NOAA SOON 
interval (1977-1998), Rome sunspot areas tended to be somewhat larger than areas measured by SOON 
(average deviation equal to 112.2 millionths of the visible solar hemisphere, with a standard deviation of 
94.2 millionths of the visible solar hemisphere). 
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Elapsed time in years from E (Rmin) 

Figure 3. Comparison of cycle 23 against mean of cycles 12-22 for R (lower panel), 
A (middle panel), and AIR (upper panel). 
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Figure 4. Comparison of s x and s R (lower panel) A x and A R (lower- middle panel), A x /s x 
and A r /s r (upper-middle panel), and n (upper panel) for 1958-1998. 


Concerning the ratios, ratios for Rome Observatory tended to be below those using RGO areas and 
Zurich sunspot numbers, but above those using SOON areas and SIDC sunspot numbers. For 1958-1976, 
the ratios based on RGO and Zurich data, on average, are ~1.5 units higher than the ratio based on Rome 
Observatory data. For 1981-1998, the ratios based on SOON and SIDC data, on average, are -2.6 units 
lower than the ratio based on Rome Observatory data. (For the brief interval of 1977-1980, when the tran- 
sitions occurred, the ratio based on Rome Observatory data, on average, was =1 .9 units higher.) 

Figure 5 displays scatterplots of area versus sunspot number, where the two plots on the left refer 
to Rome Observatory data for 1958-1976 (bottom-left panel) and 1977-1998 (top-left panel) and the two 
plots on the right refer to RGO and Zurich data (bottom-right panel) and SOON, Zurich and SIDC data 
(top-right panel). All scatterplots have linear correlation coefficients r > 0.97 and are statistically signifi- 
cant. The importance of these plots is that the Rome Observatory data remain consistent for 1958-1976 
and 1977-1998, while a sharp change occurs in the relation after the end of the RGO record. For example, 
given a sunspot number of 100 implies a sunspot area of 1,388.8 ± 134.4 or 1,363.8 ± 156.6 millionths 
of the visible solar hemisphere for 1958-1976 and 1977-1998, respectively, using Rome Observatory, 
and a sunspot area of 1,578.3 ± 120.8 or 1,081.0 ± 117.7 millionths of the visible solar hemisphere for 
1958-1976 and 1977-1998, respectively, using the RGO, Zurich, and SIDC data. 
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Figure 5. Scatterplots of A R versus s R (lower-left panel) and A x versus s x (lower-right panel) 
for 1958-1976. Scatterplots of A R versus s R (upper-left panel), and A x versus s x 
(upper-right panel) for 1977-1998. 


10 


3. CONCLUSION 


The preceding shows clearly that the relation between yearly averages of sunspot area A and sun- 
spot number R is not simply related, as has often been stated, as A= 1 6 JR. For the entire interval of 1875- 
2004, a better description is A = (13.1 ± 3.5) R, where 13.1 is the mean value of the ratio AIR and 3.5 is the 
standard deviation of the ratio. During the RGO timeframe (1875-1976), the relation is better expressed 
as A = (14.1 ± 3.2) R, while during the most recent USAF/NOAA SOON timeframe (1977-2004), the rela- 
tion is better expressed as A = (9.8 ± 2.1) R. Certainly, A versus R is highly linearly correlated (r >0.97) for 
both timeframes, although it appears very likely that A as measured today is underestimated, and R may 
be slightly overestimated. 

Figure 6 displays the variation of the minimum and maximum values of the AIR ratio for cycles 
12-23. Minimum A/R averages =8.5 ± 2.9 for cycles 12-23, 8.9 ± 3.2 for cycles 12-20 (the RGO/Zur- 
ich timeframe) and only 7.3 ±1.4 for cycles 21-23 (the SOON/Zurich/SIDC timeframe), where the first 
number in each interval is the mean value and the second number is the standard deviation. Maximum 
A/R averages, respectively, =17.0 ± 2.7, 18.4 ±1.1 and only 12.8 ± 0.9 for the same timeframes. Thus, it 
is anticipated that cycle 24 will have a minimum A/R of =7.3 ±1.4 that will occur probably in 2006 (half 
of cycles 12-23 had minimum A/R in the year of sunspot minimum and all but cycle 17 had minimum A/R 
within 1 yr either side of sunspot minimum). For the year 2004, A/R averaged 12.1 and for the yr 2005 
(through November), it has averaged 13.3 ± 4.3— inferring that minimum A/R likely will occur in 2006 or 
2007 and will have a value lower than 12.1. 

Maximum A/R for cycle 24 should be =12.8 ± 0.9, with the peak value occurring anytime during 
the maximum phase of the cycle. Thus, if cycle 24 should turn out to be =140 in size, 1:1 17 one anticipates 
a maximum yearly average for sunspot area of =1,792 ± 126 millionths of the visible solar hemisphere, 
based on USAF/NOAA SOON observations. 
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Figure 6. Variation of AIR minimum and maximum values for cycles 12-23. 
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